Synthèse sonore de clarinette avec modèle de résonateur à trous latéraux. 
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Résumé 

Les méthodes de synthèse sonore par modèles physiques 
reposent sur une description des divers éléments consti- 
tuant l'instrument. Les contraintes liées à la synthèse 
temps-réel requièrent l'emploi de géométries simples pour 
le résonateur. Cependant le corps de la clarinette est 
bien plus complexe qu'un simple tuyau parfaitement cy- 
lindrique du fait notamment de la présence des trous 
latéraux. Nous présentons ici une manière de prendre en 
compte de façon globale et simplifiée l'ouverture d'un cer- 
tain nombre de trous latéraux dans le cadre de la synthèse 
sonore en temps-réel. 
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Introduction 

Les méthodes de synthèse sonore reposant sur les modèles 
physiques nécessitent, dans un premier temps, de s'atta- 
cher à observer et comprendre quels sont les phénomènes 
les plus importants dans le mécanisme de production du 
son. Il ne s'agit pas de chercher à donner une description 
exhaustive des particularités de l'ensemble de l'instru- 
ment, mais de s'assurer que les caractéristiques des sons 
synthétisés soient proches de celles des sons de l'instru- 
ment réel. Pour la clarinette, un premier modèle sim- 
plifié de résonateur est un simple cylindre dont la lon- 
gueur correspond à la distance entre le bec et le premier 
trou ouvert. Partant de l'observation d'une différence si- 
gnificative entre les mesures d'impédances et le modèle 
cylindrique, nous nous sommes intéressés à la manière 
dont la prise en compte de la partie du corps de l'instru- 
ment placée en aval du premier trou ouvert modifie les 
propriétés acoustiques du résonateur. 

Contrairement à d'autres méthodes où chacun des trous 
latéraux est pris en compte de manière détaillée (van 
Walstijn [8] et Scavone [6]), nous nous sommes attachés à 
modéliser globalement l'effet de l'ensemble des trous ou- 
verts en nous appuyant sur les travaux de Benade [1] et 
Kergomard [2, 3]. Une formulation temporelle numérique 
du comportement du réseau de trous latéraux ouverts a 
été élaborée afin de raffiner les modèles de résonateurs 
déjà mis en œuvre dans des algorithmes de synthèse so- 
nore en temps-réel par modèle physique (Guillemain [9] ) . 

Modèle physique du résonateur 

L'étude qui suit se place dans le domaine de l'acoustique 
linéaire : les mesures d'impédance sur lesquelles sont 
visibles les différences de comportement entre le corps 
d'une clarinette et un simple tube cylindrique ont été 
présentées par Benade [1]. Elles ont été mesurées à faible 



niveau : les différences ne proviennent pas d'éventuelles 
non-linéarités. La figure 1 montre un exemple typique 
de mesures réalisées dans le cadre de cette étude et une 
impédance calculée à partir d'un modèle cylindrique. Au 
vu des dimensions et du contenu spectral du son de l'ins- 
trument, nous ne considérons que la propagation d'ondes 
planes dans la colonne d'air : la fréquence d'apparition du 
premier mode non plan pour un conduit cylindrique de 
rayon r — 7mm se trouve en effet aux alentours de 14kHz 
ce qui est bien en deçà des fréquences de résonance prin- 
cipales du résonateur. Sous ces hypothèses, la théorie de 
Kirchhoff permet de prendre en compte la viscosité et 
les effets thermiques dans le nombre d'onde k. Le rayon 
du guide d'onde étant bien plus grand que les épaisseurs 
des couches limites de pertes visco-thermiques l v et l t , 
la constante de propagation est approchée classiquement 
par : 



y i _ UJ -3/2 V^-u + (C pv — 1) Vh [Û_ 

c r V c 



(1) 



avec c = 340m/s, l v = 4 X 10" 8 m, l t = 5.6 x 10" 8 m et 
C pv = 1.4 (rapport des chaleurs spécifiques à pression et 
volume constant). 
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FlG. 1 - Impédances d'entrée réduites, mesurée sur une 
clarinette réelle pour un doigté de Fa#2 (en trait plein) 
et calculée à partir d'un modèle de résonateur cylindrique 
(en pointillés). 

Lignes de transmission 

Ces différentes hypothèses nous permettent d'adopter le 
formalisme de ligne de transmission dans le domaine 
fréquentiel (avec une dépendance temporelle en cxpjutt). 
Pour un cylindre, pression et débit acoustiques en entrée 
x = s'expriment en fonction des valeurs en sortie x = L 



pour une pulsation u) donnée : 




oùM(L)=( C0S ^ jSin ( { ^) (2) 
\j an (kL) cos (kL) J 

où Z c = pc^r 2 )^ 1 est l'impédance caractéristique de 
la colonne d'air cylindrique, avec p masse volumique de 
l'air. La dépendance de l'impédance d'entrée Z e du tuyau 
à sa charge Z s s'en déduit aisément : 

e{ ' U e (w) l+jZ s (Lo)Zë 1 ta,n(kL)' 1 ' 



son impédance itérative dès lors qu'il y a plus de deux cel- 
lules. Ceci vaut en particulier pour les basses fréquences. 
Au-dessus d'une fréquence de coupure, qui correspond, 
pour des raisons de symétrie, à la fréquence propre de 
la cellule fermée à ses deux extrémités, on a une bande 
passante. Dans celle-ci, le rayonnement devient très ef- 
ficace, et les ondes sont fortement atténuées, ce qui en- 
traîne aussi pour l'impédance d'entrée du réseau d'être 
proche de l'impédance itérative. Benade [1] avait bien 
noté ce point important, Kergomard [3] l'a interprété en 
l'associant à l'interaction extérieure des trous latéraux : 
une impédance d'entré de clarinette mesurée telle que 
représentée figure 1 montre de fait qu'il n'y a plus de 
résonances au-dessus de la coupure (environ 2kHz). 



De la cellule élémentaire. . . 

Afin de parvenir à une modélisation globale d'un réseau 
de trous latéraux ouverts, nous nous intéressons en pre- 
mier lieu au comportement d'une cellule élémentaire (Fi- 
gure 2) comportant un tronçon de longueur 2s de conduit 
principal (de rayon r) muni d'une cheminée latérale de 
rayon r t et de hauteur h c incluant les corrections de lon- 
gueurs aux deux extrémités de la longueur h [7]. 
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FiG. 2 - Schéma d'une cellule élémentaire : notations. 

La matrice de transfert entre les grandeurs en entrée et 
en sortie de la cellule correspond donc aux propagations 
sur les éléments de longueur s en amont et en aval de 
la cheminée, ainsi qu'à la dérivation au niveau du trou 
(conservation du débit et de la pression) : 
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où Y^ 1 = jZ c kh c (r/r t ) est l'admittance présentée par 
la cheminée à la colonne d'air principale en supposant 
que les dimensions du trou sont petites devant la longueur 
d'onde (nous ignorons l'impédance en série liée au champ 
antisymétrique dans la cheminée). 

... au réseau de trous latéraux ouverts 

La configuration étudiée ci-dessus ne s'applique qu'à 
un nombre limité de doigtés - ceux pour lesquels 
un seul trou est ouvert. Pour les autres, Benade 
[1] a montré expérimentalement qu'il est possible de 
considérer le réseau périodique, les trous étant identiques 
et équidistants, ce que nous discuterons plus loin. Il existe 
alors pour le réseau de trous ouverts, selon les fréquences, 
des bandes d'arrêt et des bandes passantes. Dans les 
bandes d'arrêt, les ondes du réseau sont évanescentes, 
et l'impédance d'entrée du réseau est en pratique égale à 



Sous cette condition de forte atténuation et si l'on ad- 
met que le réseau est périodique, l'impédance d'entrée 
est l'impédance itérative [2] : 



/ m 2 s 2 + ks cot (ks) „ r t 1 

Zi = jZ c tan (ks) \ / — — ; — : 71-^ ou m 



m 2 s 2 — ks tan (ks) 



r y/2sh c 
(5) 

Dans l'expression (5) apparaissent une infinité de cou- 
pures qui définissent des bandes fréquentielles aux com- 
portements variés : 

- en basses fréquences, jusqu'à la coupure proche de 
f c ~ mc/(2ir) (la plus faible pulsation solution de 
m 2 s 2 = ks tanks) de l'ordre de 1 à 2kHz, le réseau se 
comporte comme un élément de longueur équivalente 
s' = sy/l + (ms)~ 2 . L'hypothèse de troncature du 
résonateur au niveau du premier trou ouvert reste 
valable dans cette bande de fréquence et pour des 
cheminées courtes et /ou larges (ms grand). Dans le 
cas contraire, les ondes stationnaires basses fréquences 
s'étendent sur une partie du guide en aval du premier 
trou. 

- au delà, jusqu'à k — n / (2s) soit jusqu'aux environs 
de 4 à 5kHz, se trouve la première bande passante du 
réseau. 

- on observe ensuite, pour le modèle décrit par 
l'équation (5), un succession de bandes d'arrêt et de 
bandes passantes. . . 

Les considérations précédentes reposent sur l'hypothèse 
que les cellules sont identiques. En pratique, on peut 
montrer que si chaque cellule a la même fréquence de 
coupure, à condition que ses dimensions soient petites 
devant la longueur d'onde, le réseau se comporte exacte- 
ment comme un réseau périodique. Benade [1] a constaté 
que la fréquence de coupure change très peu d'un doigté 
à l'autre (hors doigtés de fourche), ce qui entraîne l'ho- 
mogénéité de timbre. Ceci semble témoigner d'un réseau 
parfaitement périodique. Pourtant, afin d'assurer une 
progression logarithmique en fréquence des notes, il faut 
écarter les trous de plus en plus quand on s'éloigne de 
l'anche, et on constate que le diamètre des trous croît 
également, assurant une fréquence de coupure à peu près 
constante. 



Le résonateur complet 

Les pics de résonance du résonateur sont modifiés par 
la présence du réseau de trous latéraux ouverts placé à 
l'extrémité d'un guide cylindrique : 

- ceux qui se trouvent dans les bandes d'arrêt du réseau 
sont légèrement déplacés vers les basses fréquences 
puisque les ondes stationnaires s'établissent sur la lon- 
gueur du cylindre augmentée d'une correction fonc- 
tion des paramètres géométriques des ouvertures et de 
l'écartcmcnt entre les trous, 

- ceux qui correspondent à la bande passante du réseau 
sont fortement atténués car une grande partie de 
l'énergie est prélevée par le réseau : les résonances ne 
peuvent plus s'établir. 

Mise en œuvre numérique 

En se focalisant sur la première bande d'arrêt et sur la 
première bande passante, il est possible d'assimiler - dans 
un domaine de fréquence de l'ordre de — 5000Hz - le 
réseau vu par le tuyau cylindrique à un système absor- 
bant les hautes fréquences, c'est-à-dire un filtre passe- 
bas pour le coefficient de réflexion présenté par le réseau 
au cylindre. C'est cet effet que nous avons cherché à re- 
transcrire de manière adaptée aux méthodes numériques 
de synthèse sonore. Une exigence qui accompagne le pro- 
cessus de synthèse sonore est le choix de la simplicité 
des modèles adoptés. Nous nous imposons l'emploi de 
filtres numériques d'ordres faibles mais dont les coeffi- 
cients sont reliés aux paramètres géométriques et acous- 
tiques du réseau. 

Il est important de noter que la principale difficulté pro- 
vient du fait que, à cause de la forme de l'équation 5, les 
résonances du réseau considéré ne sont pas des résonances 
classiques auxquelles sont associées un nombre fini de 
pôles dans le plan complexe. 

Le tracé de la réponse impulsionnclle h(t) associée à 
l'impédance réduite du réseau seul indique (Figure 3) 
la présence d'une discontinuité à l'instant initial t = 
et d'oscillations décroissantes, ce qui rend pertinent la 
décomposition de h(t) en la somme d'une impulsion de 
Dirac de hauteur ho et d'une fonction sinusoïdale amor- 
tie : 

h(t) = hoô(t) - Ae- 0t sin (uj t) (6) 
Zi(M =ho A - x2 (7) 

La contribution à l'instant initial de la partie oscillante 
semblant nulle, la hauteur ho peut être évaluée analyti- 
quement en fonction des paramètres m et s. Les coeffi- 
cients restants sont déterminés pour que le filtre présente 
le même comportement dans les deux premières bandes 
que l'impédance itérative du réseau. De plus, le réseau 
étant un élément passif (il n'est pas source d'énergie), le 
coefficient de réflexion doit présenter un module inférieur 
à l'unité. On peut montrer que ceci impose q r > 1. Si l'on 
note A les pertes visco-thermiques à la pulsation de cou- 
pure toc, les conditions précédentes mènent aux valeurs 
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FiG. 3 - Réponse impulsionnclle h(t) associée à 
l'impédance d'entrée Zi(ui) (équation (5) ) du réseau de 
trous latéraux. 

suivantes : 

h = A = f(m,s), (8) 

Qr = 1, (9) 

lu t . = cyjm 2 + A 2 . (10) 

Ce modèle analogique est ensuite discrétisé à la fréquence 
d'échantillonnage F e = 44100Hz à l'aide de schémas 
numériques centrés : 

ju> — Dz = ^(z-z- 1 ) (II) 
{juf D 2 z = F e 2 (z-2 + z- 1 ). (12) 

On obtient ainsi un filtre numérique Zi(z) introduisant 
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FiG. 4 - Coefficient de réflexion en entrée du réseau : 
modèle de réseau (en pointillés) et filtre numérique ap- 
proché (en trait plein), r = 7mm, s = fem et m = 43m" 1 . 

une atténuation des pics de résonance du tube cylin- 
drique au-dessus de la fréquence de coupure. La présence 
d'un pic d'impédance du réseau correspondant à une cou- 
pure très raide pour le coefficient de réflexion (Figure 4) 
nécessiterait, pour être prise en compte plus précisément, 
l'emploi d'un filtre numérique d'ordre très élevé. Le filtre 
retenu restitue cependant relativement bien le comporte- 
ment du réseau dans les premières bande d'arrêt et bande 
passante. 
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FiG. 5 - Impédances d'entrée du résonateur complet avec 
l'impédance itérative Zi(tu) (en pointillés) et avec la mise 
en œuvre numérique Zi(z) (en trait plein). L = 0.5m, 
r = 7mm, s = 1cm et m = 43m" 1 . 

On prend ensuite en compte le comportement acoustique 
du tuyau qui prend être modélisé comme un retard as- 
socié à un filtre numérique du 1° ordre modélisant la 
propagation avec pertes visco-thermiques [9] . Il est alors 
possible de décrire numériquement le résonateur complet 
par une équation aux différences - dont les coefficients 
sont reliés aux différents paramètres géométriques - de 
la forme : 

Pe(tn)=U e {t n ) + V (13) 

où V est fonction d'un petit nombre de valeurs passées 
et connues de la pression et du débit d'entrée. 

Optimisation 

Les valeurs des différents paramètres du modèle com- 
plet sont déterminés par comparaison à une référence. 
Nous disposons d'une série de mesures d'impédances 
d'entrée de clarinette pour une série de doigtés du re- 
gistre chalumeau, réalisées en chambre sourde au Labora- 
toire d'Acoustique de l'Université du Maine. On dispose 
ainsi des valeurs associées au mode plan [4] pour un en- 
semble de 300 fréquences entre 100 et 2500Hz. L'objectif 
est d'ajuster la longueur L et le rayon r du tuyau ainsi 
que les paramètres m et s du réseau de sorte que le filtre 
numérique Zi(z) élaboré à partir de ces valeurs coïncide 
avec la référence expérimentale Z re j pour les pulsations 
lu u considérées. Pour cela, nous avons eu recours à des 
méthodes d'optimisation globale comme les méthodes de 
recuit simulé avec une fonction coût de la forme 
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(14) 

/(w) est une fonction de pondération qui permet d'ac- 
corder une pénalité accrue aux écarts entre les pics 
d'impédance obtenus expérimentalement et ceux avec le 
modèle numérique, sanctionnant ainsi un mauvais ajuste- 
ment en fréquence et / ou en amortissement des résonances 
du système. 

L'optimisation aboutit à des valeurs de L et r relative- 
ment proches de celles mesurées sur la clarinette utilisée 



pour les mesures d'impédance, mais aussi à des résultats 
pour les paramètres du réseau de trous ouverts qui sont 
cohérents avec les dimensions et l'écartement des ouver- 
tures du corps de la clarinette. Ceci semble confirmer que 
le filtre numérique Zi retranscrit de manière efficace l'ef- 
fet du réseau de trous ouverts. 

Conclusion 

Nous avons, à l'heure actuelle, déjà pu comparer les auto- 
oscillations obtenues à partir du modèle cylindrique, du 
modèle numérique complet (cylindre et réseau) et des 
impédances mesurées (en utilisant la technique décrite 
dans [5]). D'un point de vue du timbre, les sons cor- 
respondant au modèle incluant le réseau sont beaucoup 
plus proches des sons synthétisés avec les impédances me- 
surées que ceux obtenus à partir d'un simple cylindre. 
En revanche, nous n'avons pu comparer le fonctionne- 
ment en auto-oscillation pour le filtre numérique adopté 
pour le réseau et pour un filtre présentant une cou- 
pure bien plus marquée. Il faudra de plus s'intéresser 
de manière théorique à l'influence de l'atténuation des 
pics situés au-dessus de la fréquence de coupure du 
réseau sur le fonctionnement global de l'instrument, no- 
tamment sur les régimes non-stationnaires. D'un point 
de vue de la synthèse, l'étape suivante pourra consister 
en l'élaboration d'un modèle de rayonnement corrélé au 
modèle de réseau de trous latéraux ouverts. 

Ce travail a été effectué dans le cadre de Consonnes, 
projet soutenu par l'Agence Nationale de la Recherche. 
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